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ABSTRACT 

The redshifted 21 cm transition line of hydrogen tracks the thermal evolution of the 
neutral intergalactic medium (IGM) at “cosmic dawn,” during the emergence of the 
first luminous astrophysical objects (~ 100 Myr after the Big Bang) but before these 
objects ionized the IGM 400 — 800 Myr after the Big Bang). Because X-rays, in 
particular, are likely to be the chief energy courier for heating the IGM, measurements 
of the 21 cm signature can be used to infer knowledge about the first astrophysical X- 
ray sources. Using analytic arguments and a numerical population synthesis algorithm, 
we argue that the progenitors of supermassive black holes (SMBHs) should be the 
dominant source of hard astrophysical X-rays—and thus the primary driver of IGM 
heating and the 21 cm signature—at redshifts z > 20, if (i) they grow readily from 
the remnants of Population III stars and (ii) produce X-rays in quantities comparable 
to what is observed from active galactic nuclei and high-mass X-ray binaries. We 
show that models satisfying these assumptions dominate over contributions to IGM 
heating from stellar populations, and cause the 21 cm brightness temperature to rise 
at z > 20. An absence of such a signature in the forthcoming observational data 
would imply that SMBH formation occurred later (e.g. via so-called direct collapse 
scenarios), that it was not a common occurrence in early galaxies and protogalaxies, 
or that it produced far fewer X-rays than empirical trends at lower redshifts, either 
due to intrinsic dimness (radiative inefficiency) or Gompton-thick obscuration close to 
the source. 

Key words: cosmology: theory, cosmology: dark ages, reionization, first stars, 
quasars: supermassive black holes, intergalactic medium 


1 INTRODUCTION 

Supermassive black holes ( SMBHs) reside in the hearts of 
most massive galaxies (see iKormendv &: h 3 l2013l for a re¬ 
view). Through luminous quasar episodes, they may play key 
roles in shaping their galactic and intergalactic envi ronments 
(e.g. lRicotti fc Ostrilmil[2004l : ICattaneo et al.ll200^ . Discov¬ 
eries of quasars at redshifts 2^6 — 7 have revealed that 
SMBHs with masses of several 10® Mq were alrea dy in place 
when the Universe was less thana Gyr old (e.g. |Fan et_aJj 
200l|; WiUotTWhTm rel^^ 20031: Willott et al. 2010l : 


Mortlock et al.ll 201 ll: IVenemans et al.feoiol). 

Despite the astrophysical significance of SMBHs and 
their presence throughout cosmic time, their origins are not 
well constrained by current observations an d remain a sub¬ 
ject of active investigation (see reviews bv IVolonterill20f5 : 
lHaimarJl2013l) . 


E-mail:takamitsu. tanaka@stonybrook.edu 


Broadly speaking, theoretical hypotheses of “seed” 
BHs that grow into SMBHs fall into one of two cat¬ 
egories. The first type of seeds are remnants of the 
earliest stars (Population III, or PopIII, stars), which 
form with masses 10 — 100 Mp^ at redshifts 2 > 20 
llAbel. Bryan fc Norman! |2002|; Bromm, CoDPi fc LarsonI 


200^ _ |YoshldajOmukai_&_^Ienioul^ 

Stacy. Greif fc Brommll201Cll : iGreif et al.ll20lil : lHir 


2008 


^ _.._ - - ^ . irano et al.| 

2014h . and subsequently grow thr ough a combinat i on of 


gas accretion and mergers (e.g., iHaiman fc LoebI l200ll . 
iMadau fc Reed 1200 ll : additional references in the detailed 
description in §2). The second class of seeds, called 
“direct collapse” BHs (DCBHs hereafter), forms with 
much greater initial masses, ~ lO'^ — 10® Mq at later 
times 2 < 15. They form inside gas clouds that do not 
fragment to form ordinary stars, but instead collapse 
to form a much more massive compact object (e.g . 
Bromrn_fc-Jx)glj 200 j, Knughiai>p.M,__Bullp^ 2004 


Beeelman, Volonteri fc ReesI 20061 . Lodato fc NataraianI 
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l2006h . Both families of seed models require the seed BHs to 
grow at rates comparable to the Eddington rate to explain 
the observed 10^ Mq quasar SMBHs at 2 « 6 — 7, with 
an average e-folding timescale no long er than 40 — 5 0 Myr 
between seed formation and 2 « 7 (e.g. iTanakallioIll . 

Currently, there are few observational constraints on 
the cosmic history of SMBHs at 2 > 6. The 2 ~ 6 quasar 
luminosity function and inferred underlying SMBH mass 
function leaves many degrees of freedom for theoretical ex¬ 
planations. The main empirical constraints are that mod¬ 
els should grow enough ~ 10® — 10® Mq SMBHs to ac¬ 
count for the quasar observations, and that they exceed 
neither empirical estimates of the universal SMBH mass 
density nor the unresolved cosmic X-ray back ground (e.g. 
iTanaka fc Haim^l2009l : ISalvaterra et al.ll2012li . There is no 
hrm empirical evidence that favors either kind of seed model, 
places meaningful constraints on how common the seeds 
were, or indicates whether they grew steadily (e.g., spurred 


shifts 

spurt i 


Li et al.ll2007l: Tanakall2014h or in shorter, intermittent 

^1 (IVolonteri & Reed 20051: 

Madau. Haardt & Dottil 

Volonteri. Silk & Dubudl2015 

). Even future direct ob- 


servations of quasars at 2 > 7 may not be able to distin¬ 
guish between the PopHI and DCBH seed scenarios—most 
published SMBH growth models for both scenarios are con¬ 
sistent with the existence of ~ 10® Mq SMBHs at 2 ~ 12 (a 
condition that comes about naturally if the progenitors of 
the 2 ~ 6 — 7, ~ 10® Mq quasar SMBHs grew at near the 
Eddington limit). 

In this work, we show that upcoming observa¬ 
tions of the sky-averaged redshifted 21 cm line from 
the hyperhne transition of neutral hydrogen can 
help elucidate the growth history and abundance 
of SMBHs at 2 > 20 (iRicotti. Ostriker fc GnedinI 


20051: 


_ iRipamonti. Mapelli fc Zaroubil ^ 

Mirocha. Harker fc Burns! 2013h . This iine appears 


20081 : 


absorption if the gas spin temperature is lower than that of 
the cosmic microwave background (CMB), and in emission 
otherwise; it can thus be used to map the thermal history of 
the intergalactic medium prio r to cosmic reionization (see 
iFurlanetto. Oh fc Briggsl[2006l for a review). The strength of 
the 21 cm transition is particularly sensitive to the thermal 
state of the gas, due to the coupling of the spin temperature 
with the gas temperature via c ollisions an d the Lyman-a 
background (IWouthuvsenI Il952l : iFieldlIl958l ~l. Astrophysical 
sources at “cosmic dawn,” prior to cosmic reionization 
(2 > 11), leave an imprint in this line by building up a 
Lyman-a backg round and by heating and ionizin g the 
IGM (see, e.g., 


~ound and by neatmg and ionizing tne 
Furlanett^ 20061 : Pritchard fc Foe'S I 2 OO 8 I: 


iMirocha. Harker fc Burii^ 20131 : Yaiima fc Khochfarj|201'3 L 
X-rays can contribute strongly to this signature because 
they can raise the IGM spin temperature to above the GMB 
temp erature before fully reionizing it llRicotti fc Ostrikerj 
|2004 1. The imprint should be seen in the sky-averaged 
(global) signature, because the long mean-free path of X- 
rays (with energies > 1 keV) allows them to heat the 


^ Observations do indicate that the duty cycle of quasars in¬ 
creases to > 50% at 2 > 4, compared to < 1% at 2 < 2—see, e.g., 
refs.—suggesting that SMBH growth is not very intermittent, at 
least at 4 < 2 <<. 


IGM nearly isotropically. The two strongest classes of X-ray 
sources at cosmic dawn are expected to be seed BHs ac¬ 
creting gas en route to growing into SMBHs and high-mass 
X-ray binaries (HMXBs). 

In this paper, we argue that if the SMBHs observed as 
quasars at 2 ~ 6 grew from PopIH seeds via radiatively ef¬ 
ficient gas accretion, they should leave a strong increase in 
the 21 cm brightness temperature at 2 > 20. The absence of 
this feature in future observations would imply either that 
most SMBHs formed at later epochs, that their growth was 
rare, or that they produce much less energy in X-rays rela¬ 
tive to their mass growth than the standard accretion-disc 
interpretation of luminous AGN activity. 

This paper is organized as follows. In 32 we present a 
series of analytic estimates to show that, for a wide range 
of assumptions, the X-ray output due to SMBH growth at 
2 > 6 dominates over that associated with HMXB activity. 
§ [2 presents the physical and numerical implementation of 
the PopHI and the DG models for the growth of the SMBHs, 
and the corresponding heating history of the IGM in the two 
scenarios. The computation of the 21 cm signal is presented 
in § 13.11 together with the contribution from stars. We sum¬ 
marize our main results and conclude in § [S] 


2 ANALYTIC ESTIMATES OF X-RAY 
OUTPUTS 

In this section, we compare the X-ray output from two 
classes of astrophysical sources—accreting nuclear BHs in 
the earliest haloes (galaxies), and HMXBs. 


2.1 Estimates of total emitted X-ray energy 
densities 


The comoving SMBH mass density in th e local Universe is 
Psmb h( ^ = 0) 4 X 10® M p)Mpc ~® (e.g. All er fc Richstond 


2 OO 2 I: Yu &: T remaine! 

Shankar et al.l l2004l . 


arconi et al.l 
Linking AGN luminosity 
growth 


20041 : 


func- 

tions with S MBH mass growth via Soltan’s argument 
dSoltanl I 1982 I ') , the same quantity at 2 ~ 6 is estimated 
to be psmbh( 2 : « 6) > lO"' MqMpc”®, i.e. ninety per 
cent is thought to have been accumulated at a radia¬ 
tive efficiency 77 = L/jMc?) 0.07 via gas accretion 
since 2 ~ 6 ( Sj 2 jmkar. Weinberg fc Miralda-Escudd l2009l : 


IShankar et al.l 2010ll . On average, AGN emit a fraction 
.fboi 0.05 of their total light in 2 — 10 keV X-rays (e.g. 
iHobkins. Richards fc HernauistI l2007l . and refs, therein); 
accreting stellar-mass BHs and I MBH candidates emit th e 
bulk of their light in X-rays (e.g. iFender fc Bellonil |2012| L 
a property that’s consist ent with the standard th eory of 
luminous accretion discs dShakura fc Sunvaevlll97'3 l. which 
predicts that lower-mass BHs have harder accretion spectra 
at the same Eddington fraction. 

Therefore, if the precursors to SMBHs grew via the 
same accretion mode as in the standard picture for lumi¬ 
nous AGN at lower redshifts, the total comoving X-ray en¬ 
ergy density they would have emitted prior to 2 ~ 6 can be 
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21 cm signature of SMBH growth 3 


estimated as 


eBH,x(« > 6) < 7?/bol Psmbh(2 ~ 6) c 
- 100 


(S) (^) ■=’ 


Observations show that in the absence of X-ray AGN 
activity, star-forming galaxies at low redshift produce, on 
average, 2 — 10 keV X-ray luminosities proportional to their 
star formation rate M*, with 


‘ ixM*, 


where 


£x 


> 


10^^ erg 


( 2 ) 


|Grimmj_Gil&no^_&_Sun£aCT _ _|200^ lBasinZ£che£2 


I 2 OI 3 I : Mineo. Gilfanov fc Sunvaev 20121 : iMineo et al. 2014) 


is the factor of proportionality. Most of the X-ray luminosity 
is attributable to HMXB activity, which tracks young stellar 
populations. The above relationship can also be expressed 
as a radiative efficiency (i.e. the 2 — 10 keV energy emitted 
as a fraction of the rest mass energy of stars formed): 


_ T*,X . -2 

?7*,x = ^—- = txc 
M*c2 


1.8 X 10“ 


tX,39, 


(3) 


where £x ,39 is fx in units of 10^^erg s~^Mg^yr. 

If dark matter haloes with virial temperatures > 
10"^ K (in the atomic hydrogen cooling regime) undergo rapid 
star formati on, and convert /« ~ 10% o f their baryonic mass 
to stars fe.g. iFukugita fc Peebles|l2004l i. then the total lumi¬ 
nosity density of hard X-rays emitted by star-forming galax¬ 
ies prior to 2 6 is 

e*,x(a > 6) ~ J7,,x/*^Phaio(M > M 4 ; 2 = 6) 

iio 

Mq Mpc"^ (4) 


0.8^X,39 


Above, Mi = = 10* K) is the atomic-cooling halo 

mass threshold, flb/flo « 0.15 is the ratio of the baryonic 
fraction to the matter fraction and pha.io{M > M 4 ; 2 = 6 ) « 
3x 10® M 0 Mpc“® is the 2^6 comoving mass density locked 
inside dark matter haloes above this thr eshold (computed 
using the Sheth-Tormen mass function, ISheth fc TormenI 
I 2 OO 2 I ). Note that this estimate for e*,x is dominated by stars 
in low-mass haloes just above the atomic-cooling threshold 
Mi. 

The ratio of total X-ray energy density emitted by 
SMBH progenitors to that emitted by HMXBs is then 


eBH,x(z > 6) 
e*,x(2 > 6) 


100 


/ Pbh \ 

( r* 1 

( 1* \ 

V0.07/ 

lo.osr^’®®' 

iO.W 


(5) 


It’s plausible that the earliest galaxies emitted more X- 
rays per unit mass of star formation—i.e. that they had more 
prolific HM XB activity and thus a systematically higher 
value of lyi. iBasu-Zvch et al] (l2013l ) concluded that this 
quantity is consistent with rising as oc (1 + 2 ) out to 2 « 4. 
Using the unresolved cosmi c X-ray background as a con¬ 
straint, [Djikstra^^er^^ ([20l2) ruled out £x(z) that rises more 
steeply than a power-law oc (1 -f 2 ) f or £x. 39 (z = 0) = 3 . 
Several studies jFragos et al.||2 013al 3 iHummel et al. l2014l : 
IRvu. Tanaka fc Perna. Rll2015l i' have suggested that fx may 
be a factor ~ 1 —100 higher than what is observed in lower -2 
galaxies. Even accounting for the possibility of higher £x for 
the earliest stellar populations, equation m suggests that 


eBH,x( 2 : > 6) S> t*,x(z > 6)—at 2 > 6, and that total X-ray 
production from SMBH growth should have dominated over 
that from stellar populations. 

We should also note that the above estimates for 
the SMBH mass density only account for the massive 
nuclear BHs identified in low-redshift observations. 
For example, it is possible that intermediate-mass 
(< 10^M 0 ) reside in the out s kirts of massive galax¬ 

ies (e.g. Islam. Tavlor fc SilkI 200^: O’Leary fc LoebI 


2009|;__|Micjc ^jIol_!ey-Bockelmann fc SigurdssorJ 201 ll : 


Rashkov fc Madail 2014h or in the nuclei of dwarf galaxies 


(e.g. 

|2004 


Ho et al.l I 1997 I: lizotov fc Thuaiil I 2 OO 8 I : iBarth et al.l 


Reines. Greene fc Gehal 2013lj . The expressions above 


do not account for the 2 > 6 X-ray emission from the 
growth of such BHs. 

Note that the above estimate applies to both PopHI 
and DGBH seed scenarios, as long as they arrive at the 2 « 
6 comoving SMBH mass density of ~ 10'’^ Mq Mpc“® via 
luminous gas accretion. 


2.2 Estimates of X-ray luminosity densities 

In a similar vein, we can also estimate the luminosity densi¬ 
ties (as opposed to the cumulative emitted energy densities) 
at a given redshift as a function of unknown physical pa¬ 
rameters. This quantity is more relevant for predicting the 
X-ray background as a function of redshift. 

Suppose that central BHs in galaxies possess an aver¬ 
age fraction /bh of the dark matter halo mass, and that at 
any given time they shine at an average fraction /Edd (al¬ 
lowing that a fraction of BHs are inactive, and that not all 
haloes/galaxies host a central BH). Then the 2 — 10 keV 
luminosity density of active BHs can be written as 

ZbH,X = 1.3 X 10®®/bol/BH/Edd 

X pha.io{M > Mhost) erg s“^Mq\ (6) 

where 1.3 x 10®®erg s“^Mg^ is the ratio of the Eddington 
luminosity of an object to its mass, Mhost is the minimum 
characteristic mass for a dark matter halo to host a nuclear 
BH (i.e. the halo mass scale above which the fraction of 
haloes that host a BH is close to one), and Phaio(M > Mhost) 
is the universal density of mass locked in DM haloes above 
that mass. 

The average Eddington ratio /sdd (over all haloes with 
M > Mhost, including those that do not contain an ac¬ 
creting nuclear BH) may be much higher at 2 > 6 than 
in the local Universe, where it is /Edd(2 = 0) <C 1. 
Clustering data suggest that the duty cycle of quasars in- 
creases toward high redshift, reaching ~_0.5 at 2 « 5 (e.g. 
IShankar. Weinberg fc Miralda-Escudd l2009l j. compared to 
■C 1 at 2 < 2. Such prolific activity of high-redshift SMBHs 
could be explained if th eir growth is triggered by th e major 
merger rate of galaxies et al.llioO^ : lTanakall2014l ~). which 
scales as ~ (1 -I- 2 )®^®. 

Turning to the X-ray luminosity from HMXBs, we can 
write equation @ as a luminosity density that scales with 
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4 T. L. Tanaka, R. M. O’Leary, and R. Perna 


the universal star formation rate density, 


l*,yi = O Phalo(M > M 4 ) 




C('2^)Phalo(M > M4) 

5/2 


5,7xl0»fx,(i)(l^) CW 

X Phaio(M > M4) erg s“^Mq\ 


(7) 


Above, (^{z) = —dlnphaio(Af > MT)IAz is a dimensionless 
quantity that is of order unity at redshift values of interest 
in this paper. 

The ratio of the X-ray luminosity density produced by 
accreting BH seeds to that produced by HMXBs is: 


^BH,X 

^*,x 


5 X 10 -^x.sg/Edd 




X 


V M 4 / 


l-b2 

21 


-5/2 


-1 


( 8 ) 


Above, we have used the fact that the DM halo mass func¬ 
tion dn/dMhaio roughly scales as at masses below the 

exponential cutoff of the mass function. Once again, we ar¬ 
rive at the conclusion that X-ray emission from growing seed 
BHs very plausibly dominated over the X-ray emission from 
HMXBs, unless growing nuclear BHs were rare compared 
to star-forming galaxies, or these BHs emitted much less en¬ 
ergy in X-rays per unit mass growth than their present-day 
counterparts. 

It is important to note that the above estimates hold 
even if only a small fraction of the first PopIII-forming 
minihaloes formed seed BHs with 100 Mq. This is be¬ 
cause such haloes merge rapidly, and the fraction of haloes 
containing a seed BH rapidly approaches unity (see e.g. 
iTanaka fc Haim^ 120091 '). For example, a typical halo with 
a mass > 10^^ Mq at z « 6 will have had hundreds of thou¬ 
sands to millions of progenitor haloes that formed PopIII 
stars. Below, we present results from one model that as¬ 
sumes that a X-ray luminous PopIII seed BH forms in all 
haloes that reach a virial temperature of 2000 K, and an¬ 
other model that assumes that they form in only 1% of 
such haloes. Both models arrive at approximately the same 
SMBH population by z « 6 ( cf. iTaiiaka fc Haim^ I 2 OO 9 I : 
iTanaka. Perna fc Haim^l2012i ~l. 


3 THEORETICAL MODEL 


lowed to form when the halo reaches a virial temperature 
of 2000 K, which is the approximate threshold for PopIII 
star formation. In models with DCBH seeds, we plant cen¬ 
tral BHs with masses of lO'^-lO® Mq in a small fraction of 
haloes above the atomic-cooling threshold. 

(hi) We allow BHs to grow via gas accretion, again fol¬ 
lowing semi-analytic prescriptions following plausible BH- 
halo scaling relations (motivated by empirical BH-galaxy 
scaling relations), while also satisfying observational con¬ 
straints such as the universal SMBH mass density. In our 
fiducial prescription, BHs a ccrete to approa ch a relation 
MBH(A/ha,io, 2 ) motivated bv iFerraresd (l2002l 'l . at accretion 
rates capped at the Eddingtom limit. While we treat the ac¬ 
cretion as being continuous, quantitatively similar SMBH 
populations emerge in models where BH growth occurs 
sporadica lly, e.g. triggered by major mergers of the host 
DM halo jTanaka. Perna fc Haimaril2012l : lTanakal2014l '). We 
model the X-ray emission of BHs using standard accretion 
flow theory, and assume 5 per cent of the emitted energy 
is reprocessed into a power -law corona with E > 1 keV via 
Compton upscattering (e.g. iHookins. R.ichards fc HernauistI 
l2007l : [Done et ai1l2012l ~). In all models considered in this pa¬ 
per, BH growth and HMXB activity occur only when the 
host halo is atomic-cooling (Tvir > 10“^ K). 

(iv) We allow the central BHs to merge after their 
host haloes merge (as long as the halo merger timescale 
does not exceed the Hubble time, in which case the 
haloes are assumed to become satellites), and deter¬ 
mine semi-analytically whether the merged products are 
ejected or retained after undergoing gravitational recoil (see 
ITanaka fc Haimanll2009l . for details). 

(v) Concurrently with the BH growth, we follow the uni¬ 
versal star-formation rate and the corresponding X-ray emis¬ 
sion from young stellar populations, by assuming that the 
star formation rate scales with the increase in the baryonic 
mass content of the halo (see the previous and following 
sections). 

(vi) At each timestep in the merger tree code, we track 
the cumulative X-ray background for photons with suffi¬ 
ciently long mean-free-paths (> 1 keV) to form an isotropic 
background. We account for absorption through a neutral 
IGM, as well as redshifting. We use the time-dependent X- 
ray background to compute the evolution of the mean IGM 
kinetic temperature riGM( 2 ), which in turn is used to give 
the global (sky-averaged ) 21 cm spin temperature. We follow 
IPrifAard fc Furlanettol ll2006tl and iFurlanetto fc PritchardI 
( 200fl j in evaluating the Lya background and UV emission 
from stellar populations. 


We summarize our computational scheme as follows. 


(i) Using a Monte Carlo merger tree code, we simulate 
the assembly history of dark matter haloes from 2 : ~ 50 
to 2 = 6. The algorithm ll Zhang. FakhouH_j^^Ia 200^1 and 
halo sample, which reproduce the Isheth fc TormenI (I 2 OO 2 I I 
mass func tion with high fidel ity, are the same as those de¬ 
scribed in iTanaka fc Li (l2014l l. We also account for the sup¬ 
pression of PopIII seed formation at high redshifts due to 
supersonic coherent motions of baryons against dark matter 
llTanaka. Li fc Haimanll2013l . and refs, therein). 

(ii) We plant seed BHs in DM haloes that satisfy specific 
analytic criteria, motivated by the physical models consid¬ 
ered. In models with PopIII seed BHs, the seeds are al- 


Most of the key model ingredients described above, such 
as the X-ray output from the first galaxies and the growth 
of seed BHs, contain large theoretical uncertainties (and are 
active topics of research in their own right). The goal of this 
work is not to make specific predictions for every combina¬ 
tion of plausible theoretical models, but rather to illustrate 
the fact that different families of models can correspond to 
broadly quantitatively distinct predictions for the global 21 
cm signature. Whereas previous analyses assumed that the 
growth of the X-ray background scales directly with the 
growth in the overall mass density locked inside collapsed 
dark matter haloes, here we account for relevant effects in 
the hierarchical SMBH-halo assembly process, such as in- 
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dividually limiting BH growth to the Eddington limit, oc¬ 
cupation fractions and BH-halo scaling relations, as well as 
BH mergers and the associated gravitational recoil effect. 


of X-rays are young stellar populations and (mini-) quasars, 
i.e. SMBHs (or their progenitors) that are undergoing lumi¬ 
nous accretion. Below, we describe our theoretical treatment 
of each of type of X-ray source. 


3.1 Global signatures of X-ray emission on the 
21cm transition line 


Prior to reionization, the thermal history of the IGM is most 
directly observable using the highly redshifted 21 cm line of 
neutral hydrogen. The line is globally observable in contrast 
to the temperature of the CMB with a brightness tempera¬ 
ture of 

r.”- = 4L5i„ (l - 21^) 


where xh is the neutral hydrogen fraction, Tcmb(«) is the 
CMB temperature at redshift «, and r 2 i is the spin temp er- 
ature of the neutral hydrogen (see, e.g.. lFurlanettoir200fiL for 
a more thorough discussion). The 21cm spin temperature 
of HI couples to the thermal temperature of t he gas mainly 
through the scat tering of Lyman-Q photons JWouthuvsenl 
I 1952 I : iFheldlllQS^ from the first stars, with 


TZ 




l+Xa-\-Xc 


( 10 ) 


where Xa is the coupling coefficient of the 21 cm line to 
Lyman-Q radiation and Xc is the collisional coupling coeffi¬ 
cient, whose contribution is negligible at 2 < 40. Once the 
star formation rate grows to be > lO“®M 0 yr“^, there are 
a sufficient number of Lyman-Q photons to fully couple the 
spin temper ature to the thermal tempe rature of the IGM, 
T 21 « Tigm llMcQuinn k, O’Learvll20l3 ~). 

We compute t he coefficient Xa following 
IPritchard fc Furlanettol (l2006l L summing the flux across 
the Lyman series weighted by “recycling coefficients.” The 
Lyman-Q background is computed by assuming that 9690 
photons are produced, per baryon, between Lyman-Q and 
the Lyman limit, with a power-la w spectrum oc 
(the PopII star spectral model of iBarkana fc LoebI I2005I) . 
We also calculate the dimensionless factor Sa that goes 
into computing the coeffi cient Xg, following th e wing 
approximation results of iFurlanetto fc PritchardI (l2006l : 
instead of assuming Sa = 1 ). 

Generally, the kinetic temperature of the IGM evolves 
according to 


du d f '^net 


( 11 ) 


where u, p and p are the mean specific internal energy, pres¬ 
sure and density of the IGM. The first term on the right- 
hand side, involving the derivative of the IGM density, is 
dominated by adiabatic cooling due to cosmic expansion. 
The second term includes line and continuum cooling, as 
well as Gompton heating/cooling. It also includes radiative 
heating by X-rays, which is the focus of this work. In this 
work, we calculate radiative transfer and chemistry for the 
species H, H~^, He, He'*' and He ^~*~. We refer the reader to 
iTanaka. Perna fc Haimaiil ||2012| ~) for further details. 

In a neutral IGM, photons with energies above « 1 keV 
have mean free paths longer than the Hubble horizon, and 
will build up an X-ray background that isotropically heats 
the IGM. We assume that the primary astrophysical sources 


3.2 X-ray Sources 

3.2.1 Stellar populations (HMXBs) 

Following iFurlanett^ (l2006fl . we assume that the global 
> 1 keV X-ray emission from high-redshift galaxies is pro¬ 
portional to the star formation rate (as is observed for local 
star-forming galaxies; see §2), and that as DM haloes grow 
in mass, they convert a fraction /* ~ 0.1 of their baryons 
into stars. 0 

Following eq. [2l the luminosity density in 2 — 10 keV 
X-rays produced by star formation in the early Universe at 
any redshift can be estimated as 

£2-10 keV = = I'x.*/*^Phalo(M' > M 4 ), (12) 


where Phaio(AL > M4) is the density of matter locked inside 
haloes above the atomic-cooling threshold. We assume that 
the spectral energy distribution from young stellar popula¬ 
tions takes the form of a power-law in the relevant energy 
range (i.e. har d enough to form a background), with index 
F = 1.8 (e.g. ISwartz et al.|[2003) : equation ( 1121 ) gives the 
normalization to this power-law. 

As stated in [S] we only allow HMXBs and seed BH 
growth in haloes with Tvir > 10^ K. This is a conservative 
prescription, based on the simple fact that ionization result¬ 
ing from mini-quasar activity can unb ind the gas from haloes 
with shallower potentials (see also Alyarez, Wise &: Abell 
I 2 OO 9 I : iMilosavlievic et al.ll2009l : IPark fc Ricottill 2012 l . for ad- 
ditional possible consequences of radiative feedback). 

In principle, haloes in the range 2000 K < Tvir < 10^ K 
could also produce X-rays via nuclear BH growth or HMXB 
activity. Such a scenario would result in additional heating 
as early as 2 « 30, with a corresponding turn in the 21 cm 
signature. 

It is uncertain whether the quantity which has 

been measured to be > 10 ^^ergs~^ Mlll^yr in star-forming 
galaxies in the local Universe (e. g. iGlover fc Brandi 120031 : 
iMineo. G ilfanov fc Sun vaevI 201^. evolves with redshift. 
Following IPiikstra et al.l 1 20 ldi ~ we parameterize this un¬ 
certainty by adopting a functional form 


tx,*(z} = £x,*(z = 0) X (1 -I- zf. 


(13) 


We adopt lx,*{z = 0) = 3 x 10®®ergs“^MQ^yr, and consider 
two extreme cases for the parameter fo: (i) a conservative case 
b = 0 (i.e. no redshift evolution), and (ii) a more prolific case 
b = 1, which is consist ent with observations out to 2 « 4 
dBasu-Zvch et al1l2013l L but lies near the 1 — a exclusion 
limit for higher redshifts l|Diikstra et al.l[2012l L 


^ Although we consider HMXBs as the chief emitter of X-rays 
from stellar populations, other sources, such as supernova rem¬ 
nants, could also contribute. However, such sources can be incor¬ 
porated into the quantity £x (X-ray emission per star formation 
rate), as long as they are associated with young stars. 
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3.2.2 Seed BHs 


We consider both PopIII and DCBH seed models. 

For the PopIII case, we assume that seed BHs can form 
once a dark matter halo reaches a virial temperature of 
2000 K, which is roughly the threshold required for molec¬ 
ular hydrogen to form, leading to PopIII star formation. 
Once a DM halo reaches this virial temperature, a seed 
BH forms in the mass range 10 Mq < M < SOOMq with 
a mass function and power-law slope —1.35, excluding the 
pair-instability window of 140 — 260 M 0 . While the mass 
distrib ution of PopIII sta r s rem ains a topic of active re¬ 
search jHirano et al.l |2014| . 120151 ') , we believe our choice of 
a steep power-law is conservative as only a small fraction of 
haloes form seed BHs with masses > 100 Mq. We consider 
a case where a PopIII seed BH forms in each halo reaching 
Tvir = 2000 K, and another in which a seed BH forms in 
only 1 per cent of such haloes. 

For the DCBH case, we assume that a seed BH with a 
mass of 10 ® M© forms in a small fraction of haloes reaching 
the atomic-cooling threshold of Tvir = lO"' K. 

This corresponds to a scenario in which pristine atomic- 
cooling haloes are irradiated by UV light from nearby star¬ 
forming galaxies, thereby suppressing ordinary PopIII star 
formation via photodissociation of H 2 molecules. Under such 
conditions, the gas can only cool via atomic cooling transi¬ 
tions, and collapses isothermally at a temperature ~ K, 
forming a so-called “supermassive” 10® Mq) star that 
leaves behind a BH of similar mass via general-relativistic in¬ 
stability. While the exact conditions necessary for this mode 
of BH formation, and how often they occur in nature, remain 
uncertain, recent theoretical estimates suggest that DCBH 
formation may be able to explain the abundance of the most 
luminous 2 ~ 6 quasar SMBHs (comoving number density 
of 1 Gpc~® at 2 ~ 6 ) but not the general SMBH popu¬ 
lation llDiikstra. Ferrara fc Mesingeil[2014l : lLatif et al.ll2015l : 


llnavoshi fc Tanakal 2014 1. 


We consider a somewhat generous DCBH-seeding pre¬ 
scription, in which 0.01 per cent ( 10 “'^) of all haloes reach¬ 
ing the atomic-cooling threshold within a redshift range 
20 > 2 > 10 form a DCBH of mass 10® Mq. This seed¬ 
ing prescription is sufficient to ensure that most galaxy-class 
haloes host a SMBH by 2 « 6 and at later epochs. 


3.3 Limitations 

In this work, we focus on the initial turnover feature of the 
global 21 cm signature. We do not follow the signature past 
2 « 12 , for the following two reasons. 

First, our computational algorithm tracks the mean 
thermal evolution of the IGM due to an isotropic X-ray 
background. Because the code does not track the spatial 
distribution of the dark matter haloes, we cannot follow the 
localized, expanding ionized regions that form around indi¬ 
vidual star-forming galaxies and quasars due to their UV 
emission (whereas > 1 keV X-rays propagate for ^ 1 Gpc 
before being absorbed). The strength of the 21 cm signa¬ 
ture depends on the column density of neutral hydrogen in 
the IGM; because our computational method cannot reli¬ 
ably track the latter quantity in an anisotropically ionized 
Universe, we refrain from modeling the signature into the 
epoch of UV-driven reionization, i.e. 2 < 12. 


Second, it is plausible that the heating of the 
IGM affects star for mation and (mini-) quasar activity 
in low-mass galaxies lIR lpamoii ti. Maoelli fc Zaroubil 120081 : 

iTanaka. Perna Sz HaimaiT 2012l ~). Once the IGM tempera¬ 

ture rises, gas infall into galaxies that fall be low the Jeans 
collapse scale (or the fi ltering mass scale; e.g. lGnedirJ|2000l : 
iNaoz fc Barkan3l200if ) may be suppressed, quenching star 
formation and (mini-) quasar activity in these sites. As a ref¬ 
erence, the Jeans mass scale becomes larger than the atomic- 
cooling halo virial mass threshold at Tigm 10® K. Because 
there are numerous uncertainties regarding the details of 
this global, thermally driven feedback, we do not attempt 
to model the regime where it becomes important. 

In other words, the initial rise in the global 21 cm signa¬ 
ture at 20 < 2 < 30 depends only on the heating rate of the 
IGM during this epoch, whereas the signature at later times 
depends on the progress of reionization, as well as how ef¬ 
fectively gas from a heated IGM can collapse into low-mass 
haloes. We concentrate on the former signature. 


4 RESULTS 

4.1 IGM kinetic temperature 

In Figure [T1 we show the redshift evolution of the mean 
IGM kinetic temperature. We will describe the curves in 
order, roughly from bottom to top, as follows. The solid 
black fine at the very bottom shows, for reference, the null 
case with adiabatic cooling due to cosmic expansion and no 
heat sources. 

With the exception of the adiabatic case, all of the mod¬ 
els include X-ray heating from HMXBs. In models with thin 
fines, it’s assumed that the factor £x,* (which relates X-ray 
production to star formation rate) is the same as measured 
in star-forming galaxies the local-Universe; in those with 
thick fines, this quantity increases toward high redshift as 
^X,t oc (1 -I- 2 ). 

The thin, magenta line and the thick, red fine (both dot¬ 
ted) show cases where HMXBs are the only X-ray sources 
heating the IGM. The former, which assumes that X-ray 
production per star formation rate (the quantity £x,») does 
not evolve with redshift, results in the IGM kinetic tempera¬ 
ture increasing at 2 ~ 20. For the latter case, which assumes 
that £x,* oc (1 -I- 2 ), Tigm increases at 2 ~ 25. 

The thick, green curve and the thin, dark green 
curve (both long-dashed) show models where DGBHs 
with M = 10®Mq form in 0.01% of atomic-cooling 

haloes at 2 < 15, and grow to match the pre¬ 

scribed relationship MBH(A4haio, 2 ), with accretion capped 
at the Eddington rate. Note that this seeding fraction 
of 0 . 01 % is much larger than found in recent theoret¬ 
ical work, which suggest that DGBHs may have chal¬ 
lenges in matching the observed abundance of h igh-redshift 
quasa rs jDijkstra, Ferrara fc Mesingeil l2014l : iLatif et al.l 
I 2 OI 5 I : IlnaTOshi fc Tanakal l2014l 'l. These heating curves are 
the same as HMXB-only cases described in the previous 
paragraph, except for the additional heating at late times 
(2 < 15) due to the emergence of DGBH seeds. 

Finally, the blue and cyan lines (both thick and short- 
dashed) show cases in which a PopIII seed BH forms in 
haloes reaching virial temperatures Tvir = 2000 K, and 
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Figure 1. The IGM temperature (as contrasted to pure adia¬ 
batic cooling) resulting from different sources of heating: direct 
collapse BHs (green and dark green lines) and BHs from PopIII 
seeds (blue and cyan lines), both including the contribution from 
HMXBs. The effect of HMXBs alone is also displayed separately 
(magenta and red lines). Two prescriptions (compatible with in¬ 
dependent constraints) are considered for each heating source. 
Heating from PopIII BHs dominates over that of HMXBs alone, 
while the contribution from DCBHs is comparable or subdomi¬ 
nant with respect to that of HMXBs, until z ^ 15. For reference, 
we have also plotted what the IGM temperature would be with 
adiabatic cooling alone (solid black curve), as well as the CMB 
temperature (which determines when the 21 cm line goes into 
emission; yellow-and-black dashed curve). 

promptly begins to accrete gas as the host halo also grows. 
The cyan line shows the case where such a growing seed 
forms in all Tvir > 2000 K haloes, and the blue line shows 
the case where a seed forms in 1 per cent of such haloes. 
Both cases assume an HMXB contribution £x,* oc (1 + 2 ), 
but it is clear that the mini-quasar activity of PopIII BHs 
dominate. 

The PopIII seed models heat the IGM much earlier than 
models that include HMXB heating only (which, again, use 
the growth in ph aiof-M > M 4 ) as a proxy for star-formation; 
lFurlanettal200d ) . Again, the former assumes that seed BHs 
form in > 1% of Tvir ~ 2000 K haloes, and grow to match 
extrapolated BH-halo scaling relations in all haloes with 
Tvir > 10^ K. 

It is important to point out that even if PopIII rem¬ 
nants produce copious X-rays, they may not contribute sig¬ 
nificantly to the present-day X-ray background. This is be¬ 
cause their spectra are expected to be relatively soft, as well 
as sharply peaked in intrinsic X-rays. For example, a black- 
body, Eddington-rate accretion disc around a 10® Mq BH 
has a spectrum that peaks at « 0.3 keV and drops off by 
orders of magnitude above » 1 keV. A graybody disc llBlaed 
|2004 iTanaka fc Menoril l2010li can be harder and peak at 
« 1 keV. In either case, the bulk of the emission would 


[t] 


redshift 

15 20 25 30 



100 90 80 70 60 50 
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Figure 2. The brightness temperature of the 21 cm radiation (top 
panel), and its spectrum (bottom panel) for a medium heated by 
HMXBs alone, or with the inclusion of accreting BHs either in the 
DCBH scenario or in the PopIII formation scenario. Line types 
and colors refer to the same models as in Fig. [T] 

be redshifted to energies well below the « 0.1 keV limit 
of soft X-ray observatories. Therefore, it is plausible that 
the vast majority of the X-rays emitted by miniquasars at 
2 : > 20 is redshifted into the FUV, and be minimally re¬ 
flected in the present-day X-ray background. Note the con¬ 
trast to present-day AGN, whose intrinsic spectra peak in 
the UV and whose X-ray emi ssion consists of a power-law 
tail (cf. ISalvaterra et al.ll2012l . who consider limits on early 
SMBH growth by assuming power-law X-ray spectra). 


4.2 21 cm brightness temperature 

After the first stars form, the thermal history of the IGM 
is directly imprinted in the 21cm brightness temperature. 
One of the most important characteristics of the global, sky- 
averaged 21 cm signal is when it reaches its minimum. This 
turning point occurs when the IGM transition s from pri¬ 
maril y adiabatic cooling to the e p och of heating jFurlarietto 
2001; IPritchard fc LoeH I 2 OO 8 I : iMirocha. Harker fc Burns 
2013I) . which is sensitive to the accretion history of the 
first black holes. A number of experiments have been pro¬ 
posed, or are currently underway, to detect this signal— 
e.g. t he Dark Ages Radio Explore]0 (DARE; iBurns et al.l 
|2012|). the Large-Aperture Experiment to Detect the Dark 
Ageu (LEDA) , the Experiment to Detect the Global EoR 
Ste p (EDGES: iBowmam Rogers fc Hewit^l2008l) . and SGI- 
HI (IVovtek et al.ll2014l) . 


^ http://lunar.Colorado.edu/dare/ 

^ https://www.cfa.harvard.edu/LEDA/ 
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The upper panel of Figure [2] shows the 21 cm bright¬ 
ness temperature as a function of redshift for the 

same models shown in Figure [T] These models use the same 
line styles and follow the same order as discussed in § 14.11 
Because PopIII BH models rapidly accrete gas, they begin 
heating the IGM as soon as they form. The total amount of 
X-ray heating dominates over stellar sources (HMXBs and 
supernova remnants). As a result, our PopIII BH models 
increase in brightness temperature beginning at z < 25, at 
typically « 20 mK above the adiabatic case. The detection 
of such a signature would strongly suggest either (i) stellar 
populations from the earliest galaxies produced far more X- 
rays per unit mass formed in stars than what is observed 
from nearby star-forming galaxies, or (ii) nuclear BHs were 
prolifically emitting X-rays as mini-AGN at these redshifts 
(much earlier than predicted by most DCBH scenarios). 

In contrast to the PopIII models, the DGBH models do 
not produce sufficient X-rays to signihcantly heat the IGM 
before z > 20. HMXBs produced in ordinary stellar evo¬ 
lution dominate the IGM heating while the 21 cm signal 
is in absorption. Because the gas has more time to adia- 
batically cool, the resulting signal is intrinsically stronger 
and easier to detect with a lower fo reground sky tempera¬ 
ture dFurlanetto. Oh fc Briggj|200d f . If the 21 cm trough is 
detected at late times, either the PopIII BHs were heavily 
obscured, or the BHs that are the progenitors of SMBHs 
formed via direct collapse. 

At the frequencies required to detect the dark ages 
(50 —90 MHz) the observations will be dominated by smooth 
Galactic foreground emission. The 21 cm signal contributes 
only about 1 part in lO'^ to the total signal. The 21 cm 
signal is detectable by the variations of the brightness tem¬ 
perature with frequency. The lower panel of figure [5] shows 
the derivative of the brightness temperature, dT^^^/di^, as 
a function of z, for the same models. 

An upward trend in the 21 cm brightness temperature 
prior to z > 20 strongly implies rampant PopIII seed growth, 
or a much higher X-ray production from young stellar popu¬ 
lations than seen in the local Universe. At these frequencies, 
planned experiments are dominated by systematic uncer¬ 
tainties. For example, DARE will have a measurement un¬ 
certainty of order 1 mK at these frequencies after 3000 h of 
integration. In addition, it is capable of determining when 
the brightness temperature reaches a minimum to 5z « 0.4. 
Typical uncertainties in the minimum temperature are ex¬ 
pected to be cm/ y21 cm ^ 15 % for niodels similar to the 

ones considered here llHarker et al.l[201^ . 

Finally, more stringent constraints on the progenitors 
of SMBHs may be made by combining the global 21 cm 
observations with other techniques. For example, observa¬ 
tions of the evolution of star formation at high redshift 
may be able to constrain the contribution of stellar sources 
to the IGM heating. Angular correlations in the 21 cm 
signal dP ritchar d fc Loe^ l2008l : iMcQuinn fc O’Learvl[201^ : 
IVisbal et al.ll2012l ') may also be able to ascertain the char¬ 
acteristic dark matter halo masses that dominate the IGM 
heating. 

Both reionization and thermal feedback on baryonic 
structure formation would suppress the 21 cm brightness 
temperature at late times. An upward trend in the bright¬ 
ness temperature at z < 15 would therefore strongly imply 
that SMBH growth began in earnest at this epoch. 


5 SUMMARY 

IGM heating by the first radiation sources leaves a marked 
imprint on the 21 cm radiation. Here, using both analyti¬ 
cal arguments as well as Monte Carlo simulations, we have 
argued that the primary driver of IGM heating are the BH 
remnants of PopIII stars, if these are the progenitors of the 
SMBHs observed in the Universe. In this scenario, heating 
by these sources would dominate over the other main con- 
stribution, i.e. HMXBs, unlike the case in which SMBHs 
grow from DCBHs, which we hnd to be subdominant with 
respect to HMXBs. 

Our work in estimating the 21 cm signature of 
the progenitors of SMBHs is distinct in three respects. 
First, the total BH mass density in our models does not 
grow in proportion to the total mass density locked in- 
side dark matte r haloes (as in most pre v ious works, e.g . 
iFurlanett^ l2006l : IPritchard fc Loe9 l2008l : iMirochal 1 201 dl l . 
but are grown individually while reflecting effects such 
as the Eddington accretion limit, mergers, and gravi¬ 
tational recoil. Second, we explicitly treat the propaga¬ 
tion and absorption of the X-rays averaged over the cos¬ 
mological volume (iMiroc h3 I 2 OI 4 I ') . Finally, our models 


are constrained to those 

that reproduce the pro 

parties 

of high redshift quasars ( 

Tanaka. Perna & Haiman 

I 2 OI 2 I: 

iTanaka. Li & Haimanll2013 

). 


The 21 cm signal, with its sensitivity to the heating 
history of the Universe, hence becomes a powerful tool to 
probe the progenitors of SMBHs, and their growth history. 
In particular, we have shown that a trend upward in the 
brightness temperature of the 21 cm radiation would indi¬ 
cate an abundant production of PopIII seed BHs, and hence 
would lend support to a scenario in which the growth of the 
SMBHs originated from those. 
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